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Abstract 

Adler and van Moerbeke [13] described a reduction of the 2D-Toda hierarchy called 
Toeplitz lattice. This hierarchy turns out to be equivalent to the one originally described 
by Ablowitz and Ladik [T] using semidiscrete zero-curvature equations. In this paper 
we obtain the original semidiscrete zero-curvature equations starting directly from the 
Toeplitz lattice and we generalize these computations to the matrix case. This gene- 
ralization lead us to the semidiscrete zero-curvature equations for the non-abelian (or 
multicomponent) version of the Ablowitz-Ladik equations [3]. In this way we extend the 
link between biorthogonal polynomials on the unit circle and the Ablowitz-Ladik hierar- 
chy to the matrix case. 
MSG classes: 05E35, 37K10. 

1 Introduction. 

The Ablowitz-Ladik hierarchy has been introduced in 1975 |1] as a spatial discretization of 
the AKNS hierarchy. As described by Suris in [2] Ablowitz and Ladik replaced the celebrated 
Zakharov-Shabat spectral problem 

< 

d^^ = AW 

with a discretized version of it; namely 

< 

dr^k = Mk^k- 

Here ^ and are two-component vectors while L,M,Lk and Mj. are 2x2 matrices. In 
particular 

' Hi 

and 

\yk z y 



We can consider, as usual, the periodic case {k G the infinite case (k E Z) or the semi- 
infinite case {k G N). The standard zero-curvature equations for the AKNS hierarchy are 
replaced with the semidiscrete zero-curvature equations 

drLk = Mk+iLk - LkMk. 

As an example, one of the most important equation of this hierarchy is the discrete complex- 
ified version of the nonlinear Schrodinger equations 

{drXk = Xk+l - 2Xk + Xk-l - XkVkiXk+l + Xk-l) 
drUk = -yk+i + 2yfc -yk~i + xkUkiuk+i + yk-i)- 

Quite recently different authors (see [13] and [15]) worked on the link between biorthogonal 
polynomials on the unit circle and the semi-infinite Ablowitz-Ladik hierarchy. This is an 
analogue of the celebrated link between Toda hierarchy and orthogonal polynomials on the real 
line. Many articles have been written about orthogonal polynomials and the Toda hierarchy 
and it is almost impossible to recall all of them. Let us just mention the seminal paper of Moser 
[H] and the monography [9]; there interested readers can find many references. The more recent 
articles [11], [12] and [13] should be also cited. In these articles the connection between the 
theory of orthogonal polynomials and the Sato's theory of infinite Grassmannians appeared 
for the first time. In particular in [13] the case of Toda and Ablowitz-Ladik hierarchies are 
treated in a similar way as reductions of 2D-Toda. In this way it's given a clear and unified 
explanation of the role played by orthogonal polynomials in the theory of integrable equations. 
As noted by the authors their approach is quite different from the original one; actually in 
their paper [13] they always speak about Toeplitz lattice and the coincidence with Ablowitz- 
Ladik hierarchy is just stated in the introduction. Nevertheless, as explained in the second 
section of this paper, it is very easy to deduce semidiscrete zero-curvature equations starting 
from Adler-van Moerbeke's equations. Having this approach in mind we addressed a question 
arising from the following facts: 

• Time evolution for orthogonal polynomials on the real line leads to the Toda hierarchy. 

• Time evolution for biorthogonal polynomials on the unit circle leads to the Ablowitz- 
Ladik hierarchy. 

• Time evolution for matrix orthogonal polynomials on the real line leads to the non- 
abelian Toda hierarchy (see for instance [17j and references therein; our approach is 
close to [20]). 

What about time evolution for matrix biorthogonal polynomials on the unit circle? 

In other words our goal was to replace the question mark in the table below with the 
corresponding hierarchy. 





Orthogonal polynomials on M 


Biorthogonal polynomials on 


scalar case 
matrix case 


Toda 
non-Abelian Toda 


Ablowitz-Ladik 

? 
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In this article we prove that the relevant hierarchy is the non-Abelian version of the 
Ablowitz-Ladik hierarchy. This hierarchy has been already studied by different authors since 
1983 (see [3], [1] and [5]) but, at our best knowledge, the connection with matrix biorthogonal 
polynomials on the unit circle was never established before. We also remark that this hierarchy 
is usually called matrix, vector or multicomponent Ablowitz-Ladik. Here we prefer to call 
it non-Abelian to stress the analogy with Toda. In our setting, instead of having one Lax 
operator of size 2 x 2, we have two Lax operators C\ and C^. of size 2n x 2n. These 
operators depend on matrices and of size n x n. Semidiscrete zero curvature 

equations are given by 

where also A^^ and Ai^. are block matrices. For instance we have, in this hierarchy, two 
versions of the non-Abelian complexified discrete nonlinear Schrodinger: 

drVl = -Vk+i + '^Vk - Vk-i + Vk+iAVk + VkAVk-i 
UrXf^ — J.f^^^ ^-f-fc "T -^fe-l -^k-ltlk-^k -^kyk-^k+l 

[drvi = -yi+1 + '^yi - vLi + yi-i^lvi + vi^lvi+i 

(see [H] for a review about these equations). 
Sections are organized as follows: 

• In the second section we recall some basic facts about 2D-Toda and the connection with 
biorthogonal polynomials. We use the approach developed in [13j . 

• Third section starts with the description of the Toeplitz lattice (see |T3| ) and shows how 
to deduce semidiscrete zero-curvature equations for the Ablowitz-Ladik hierarchy. 

• In the fourth section we extend the Toeplitz lattice to the case of block Toeplitz matrices. 

• Fifth section gives recursion relations for matrix biorthogonal polynomials on the unit 
circle; these formulas slightly generalize formulas contained in [i9\ and ^22j for matrix 
orthogonal polynomials on the unit circle. 

• In the sixth section we derive block semidiscrete zero-curvature equations defining the 
non-Abelian Ablowitz-Ladik hierarchy. As an example we write the non-abelian ana- 
logue of the discrete nonlinear Schrodinger. 



2 2D-Toda; linearization and biorthogonal polynomials. 

In this section we recall some basic facts about the 2D-Toda hierarchy as presented in [7]. 
Moreover we describe the connection with biorthogonal polynomials as originally presented 
in [13j . We are interested in the semi-infinite case; we start denoting with A the shift matrix 

A := {^i+l,j)i,j>0- 
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For the transpose we use the notation A"^ = A ^. Then we define two Lax matrices 

'Li :=A + E.<o«f^A* 

< 

M :=aSA-i + E.>oaPA^ 

where {(if \ s = 1, 2} are some diagonal matrices. 2D-Toda equations, expressed in Lax form, 
arise as compatibihty conditions for the following Zakarov-Shabat spectral problem: 







Lin 


= z-'n 




= (^?)+* 


dt^n 






= (^^)-^ 


dsr^n 


= -iL-2V. 



Here we introduced two infinite sets of times {ti,i > 0} and {si,i > 0}. We denoted with 
the upper triangular part of a matrix N (including the main diagonal) and with the lower 
triangular part (excluding the main diagonal), and are semi-infinte column vectors of 
type 

M/i(z) = (*i,o(z),*i,i(z),...)T 

mz) = m,o{^),n,i{z),...)^. 

For every k the two expressions e" 

-?(t,^)^^_^(^) and e'^^^'^'^nk 

(z ^) are polynomials in z of 

order k. Lax equations are written as 



dtM=[iL'lh,L,] ds„U=[m).,u 



1,2. 



2D-Toda equations can be linearized as explained in [7] . We start with an initial value matrix 
M(0, 0) = {Mjj(0, 0)}jj>o and we define its time evolution through the equation 

M{t;s) :=exp (^^(t. A)) M(0, 0) exp - ^(s, A"^)) . 

We assume that there exist a factorization 

M(0,0) =5i(0,0)-i52(0,0). 

Here Si is lower triangular while 5*2 is upper triangular. We assume that both Si and ^2 
have non zero elements on the main diagonal and we normalize them in such a way that every 
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element on the main diagonal of Si is equal to 1. Moreover we consider values of t and s for 
which we can write 

M{t,s) = Si{t,s)-^S2{t,s) (1) 

with Si and S2 having the same properties as above. It can be proven that such factorization 
exists if and only if all the principal minors of M{t, s) do not vanish. In particular this 
condition is satisfied when M{t, s) is the matrix of the moments of a positive-definite measure. 
Now we denote with x{z) the infinite vector x{z) '■= (1, z, z^, . . .)^. Wave vectors for 2D-Toda 
and Lax matrices are constructed in the following way. 

Theorem 2.1 ([7J). The wave vectors 

^i(z) :=eMi{t,z))Six{z) 
%{z) :=exp(-e(s,z-i))(S2-i)Tx(z-i). 

and the two Lax operators Li := SiAS^"^ and L2 := 82^^^ satisfy the 2D-Toda Zakharov- 
Shabat spectral problem. 

Proof We just sketch the proof and make reference to the article It is clear that the 
matrix M(t, s) satisfies differential equations 

dt.M = A'M 
ds-M = -MA-\ 

Then it is easy to deduce Sato's equations 

dt^Sr = -(L7)_Si 

dt„S2 = {L^,hS2 

a,„Si = (L5)_Si 

ds„S2 = -iL^) + S2. 

and Zakharov-Shabat's equations can be deduced from the expression of wave vectors in terms 
of 5i and 82- □ 

The last thing we need is the link between the factorization of M and biorthogonal poly- 
nomials. We introduce a bilinear pairing on the space of polynomials in z defining 

< z\z^ >M-= Mij. 

The following proposition is a direct consequence of ([I]). 
Proposition 2.2 ([13j). 

are biorthonormal polynomials with respect to the pairing <,>m; i-e. 
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3 From the Toeplitz lattice hierarchy to the semidiscrete zero- 
curvature equations for the Ablowitz-Ladik hierarchy. 



In this section we briefly recall the reduction from 2D-Toda to the Toeplitz lattice as described 
in . Then we show how the Ablowitz-Ladik equations are easily obtained from the Toeplitz 
lattice. Suppose that our initial value M{0, 0) is a Toeplitz matrix; i.e. we have 



M(0,0) =r(7) 















^(2) 


^(1) 


^(0) 



...\ 



■••/ 



for some formal power series ■~f{z) = T*^"^-^"- Since A = T{z ^) is an upper triangular 

Toeplitz matrix it follows easily (see for instance [23]) that 



M{t,s) = exp [C {t, A) j M(0, 0) exp - ^ (s, A"^ 
r( exp {at, z'^))jiz) exp ( - ^(s, z))) . 



This means that Toeplitz form in conserved along 2D-Toda flow, hence we are dealing with 
a reduction of it. This reduction is called Toeplitz lattice in [13]. In that article the authors 
noticed, in the introduction, that this is nothing but the Ablowitz-Ladik hierarchy. Now we 
will describe how to obtain the original formulation of the Ablowitz-Ladik equations starting 
from Adler-van Moerbeke's formulation. The key observation is that, in this case, the bilinear 
pairing < p,q >m between two arbitrary polynomials is given by 



<p,q >M= 



p{z)j{z)q* [z] 



dz 
2iTiz 



Here the symbol of integration means that we are taking the residue of the formal series 
p{z)'y(z)q* {z) and q*{z) = q{z~^). In other words q^^^ and q^^^ are nothing but orthonormal 
polynomials on the unit circle. We also define monic biorthogonal polynomials 



p 



p 

(2) 



(1) 



with h = diag(/io, /ii, /i25 • • •) some diagonal matrix. Given an arbitrary polynomial q{z) of 
degree n we define its reversed polynomial q{z) := z^q*[z) and reflection coefficients 



Xn:=p'i\Q) yn:=p'nKQ). 

We can state the standard recursion relation associated to biorthogonal polynomials on the 
unit circle. The equation below was already known to Szego ^18j. It has been used in |14j in 
relation with the theory of integrable equations. 



(2), 
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Proposition 3.1. The following recursion relation holds: 



Cr. 



(2) 



Using this recursion relation Adler and van Moerbeke in [13] wrote the pecuhar form of 
Lax operators for the ToepHtz reduction. 

Proposition 3.2 ([13j). Lax operators of Toeplitz lattice are of the following form: 



(-xiVQ 1 - xiyi 
-X2yo -X2yi 1 



h~^Lih 



-xsyo 



\ 



-xsyi 



\ 

- X2y2 
-2:32/2 1 - xsys 



/ -xoyi 
1 - xiyi 




V 



-xoy2 
-xiy2 
1 - 2:2^2 





-xoys 
-xiy3 

-X2y3 

- xsys 




From this proposition it fohows easily the following corollary. 
Corollary 3.3. Reflection coefficients satisfy the following equation 

hn+i 



hr. 



1 - Xn+iyn+\- 



Now we can state the theorem relating the Toeplitz lattice to the original form of the 
Ablowitz-Ladik hierarchy. 



Theorem 3.4. The Toeplitz lattice flow can he written in the form 



di 



for some matrices Mti,n, -^Si,n depending on {xj,yj, z}. 



^P^n\z) 



P^n\z) 



P^n\z), 

'p';l\z) 

p'nHz)^ 



(3) 
(4) 



Proof We prove it for the set of times denoted with t. 

We denote d{[z]) = diag(l, z, z^, z^, . . .). We have the identities 

^i{z)=exp{^{t,z))p^^\z) 
%iz) = h-^d{[z-^])eM-as,z'W\z) 
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giving the following time evolution for orthogonal polynomials 

dup^'\z) = -hd{[z]){V,)l^h-'d{[z-'])p''^\z). 

Here {L\)j^j^ denotes the strictly upper diagonal part of L". The formulas above are obtained 
from a straightforward computation and using the fact, proven in [T3], that 

dt^ \og{K) = {L\)nn- 

Hence we have that, for every k, dt^p^f}^ is a linear combination of {p^k \ P^^+i^ • • •} '^ith 

(2) 

coefficients in C[xj,yj]. In the same way, for every k, dt^p\ is a linear combination of 

f2) f2) f2) 

{Pfc ;P^_iiPa:_2'---} '^it^ Coefficients in C[xj,?/j]. Using the recursion relation ^ and its 
inverse 

we can obtain the desired matrices Mti,n- D 




Corollary 3.5 (Ablowitz-Ladik semidiscrete zero-curvature equations). The matrices Cn 
satisfy the following time evolution 

dt^Cn = Mt,,n+l^n " ^nMt,,n (5) 
ds.Cn = Ms„n+lCn - CnMs„n- (6) 

Proof These equations are nothing but the compatibility conditions of recursion relation ^ 
with time evolution ([3]) and □ 



Remark 3.6. Actually our Lax operator £„ is slightly different from the Lax operator L„ used 
in Ul by Ablowitz and Ladik and written in the introduction above. Nevertheless, as shown 
in [16], these two Lax operators are linked through a simple change of spectral parameter. 

Example 3.7 (The first flows; discrete nonlinear Schrodinger.). The first matrices ■M.ti,n and 
Msi,n are easily computed. We have 

a (1) (1) (1) , (1) (1) , ~(2) 

OtiPk = -^Pk - Xk+iVkPl + Pk+i = -Xk+iVkPl + Xk+ipl 



o ~(2) hk+l .(2) (1) -(2) 

OtiPl = -z-f^PlU = zykPl -zpl' 



giving immediately 



M 



-Xk+lVk Xk+l 

zyk -z 

An analogue computation for si can be easily done obtaining 



Ms^^k 



Z —Z ^Xk 

-yk+i xkyk+i 
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We can already write, with ti and si, the well known integrable discretization of nonlinear 
Schrodinger. We just need two more trivial rescaling times introduced with substitutions 

p^^^ I — > exp 



.(2) _ ...(-,o)p(2) 



and corresponding to matrices 



1 



-1 



Now we can construct the matrix 

Mr,k = Mti,k + Ms-,,k - Mto,k - Mso,k 

- \- Xk+lVk Xk+l - Z~'^Xk 

zyk - Vk+i XkVk+i + 1 - ^ 
associated to the time t = ti + si — to — sq so that semidiscrete zero-curvature equation 

dr^k = MT,k+l^k — ^kMr,k 

is equivalent to the system 

' drXk = Xk+l - 2Xk + Xk-l - XkVkixk+l + Xk-l) 

(7) 

^drVk = -Vk+i + %fc - Vk-i + XkVkiVk+i + Vk-i)- 

This is exactly the complexified version of the discrete nonlinear Schrodinger. Rescaling 
r I— > and imposing j/j, = itx^ we obtain 

- idrXk = Xk+l - 2xk + Xk-l T \\xkf{xk+i + Xk-i). (8) 

4 Toda flow for block Toeplitz matrices and the related Laix 
operators. 

Now wc generalize the Toeplitz lattice's equations to the block case. We start with a matrix- 
valued formal series 



Here every element 7^'^) is a n x n matrix. Then we define its time evolution as 

l{t, s; z) := exp ( - ^(s, z~h))-i{z) exp (^(i, zl)) . 
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where I is the n x n identity matrix. Differently from the scalar case we don't consider just 
one Toeplitz matrix but the two block Toeplitz matrices, right and left, given by 





^(-1) 


^(-2) 








^(2) 






\ ■ 










-y(2) 










^(-1) 


^(0) 



V 
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In this way we obtain the following linear time evolution for our block Toeplitz matrices (in 
the following we will omit the symbol 7): 



(9) 
(10) 



where, in this case, we have A = T^(z ^I). Then we assume that there exist two factorizations 



Si ^82 



Z2Z1 



Here are block-lower triangular while 5*2, ^1 are block-upper triangular. We assume 

that all these matrices have non degenerate blocks on the main diagonal (i.e. these blocks 
must have non zero determinants). Normalizations are chosen in such a way that every 
element on the main block-diagonal of 5i and Z2 is equal to the identity matrix I. As we 
did before wc assume that these conditions hold when every time is equal to and the we 
consider just the values of t and s for which these conditions still hold. In the matrix case we 
can define two bilinear pairings given by the following definition. 



Definition 4.1. 



<P,Q >r:= j P*i^hi^)Qi^)^ <P^Q >i-= / P{zMz)Q*iz)-, 



dz 



2mz J 2'!riz 

where P and Q are two arbitrary matrix polynomials and P*{z) := {P{z~^))^ 

Our two factorizations give exactly biorthonormal polynomials for <,>,• and <,>;. In 
the following we denote x{z) := (I, zl, z'^l, z^l, . . .)^. 
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Proposition 4.2. 

:= I Q?' I = SM^) (11) 



Qi-y := Q?' = (S^Txiz) (12) 

Q^'^'^ ■■= {q^^^' ..)=x{zfZ, (13) 

Q(^)-=(q(^> Qf)^ . . .) = X(^)^(^2-^)^ (14) 

are the biorthonormal polynomials associated to the pairing <, >; and <,>r- Hence for every 
i,j we have 

Proof We just prove, as an example, the proposition for the right polynomials. On the other 
hand the one for left polynomials is identical to the usual proof for 2D-Toda. We have 

( < Qf ^^ Qf^'' >r] = ( y: {Z^%^ < Z\ A >r (Z,),, 




z-^T^Zi = I <^ = Z2Z-^. 

(it should be noted that, in this case, subscripts of type {Zi)ij denote the block in position 
{i,j) and not the element {i,j)-) □ 

Now we can write the related Sato's equations for Si and Zj. It is convenient to introduce 
the following Lax operators. 

Definition 4.3. 

Li := SiAS^^ L2 := ^sA'^^s"^ (15) 
Ri := Z^^A-^Zi i?2 := Z^^AZ2. (16) 

Proposition 4.4. The following Sato 's equations are satisfied. 

dt^Si = -(L«)-5i = -Z^{R^)+ (17) 

dt^S2 = {L^)+S2 dt^Z2 = Z2{R'l)- (18) 
= (L^)_5i 5,„Zi = Zi(i?^)+ (19) 

ds^S2 = -{L^2) + S2 ds,,Z2 = -Z2{R'i)- (20) 

Proof Wc will just prove, as an example, the equations involving i-derivative of Z\ and Z2. 
We assume, as an ansatz, that we have 

dt^Zi = ZiA 

9t„^2 = Z2B. 



11 



for some matrices A and B. Then exploiting time evolution of we can write 

r-A-^ = dt„T^ = dtjZ2Z^') = 
ZaSZf ^ - Z2Z^^ZiAZ^'^ = Z2{B - A)Z^^ 

hence we must have {B — A)Z^^ = Z^^h~'^. We rewrite it as 

B-A = ^A-"Zi = R^. 

Moreover we have that B must be strictly block-lower triangular and A block-upper triangular. 
This is because Z\ is block-upper triangular and Z2 is block-lower triangular with constant 
entries on the main diagonal. Hence we conclude that A = — and B = {K^)-. □ 

Now it's just a matter of trivial computations to write down the corresponding Lax equations 
for Li and Rj. 



Proposition 4.5. The following Lax equations are satisfied: 

9t„Li= {Li)+,Li dt„Ri = Ri,{Ri)- 
dsn^i = {L2)-,Li ds^Ri = Ri,{R2)-\ 



(21) 
(22) 



The definition of our Lax operators will give us eigenvalue equations for suitably defined 
wave vectors. 



Definition 4.6. 

^^1(2;) :=exp(e(t,2;I))5ix(2;) 
$i(z) :=expiat,zI))[x{z)Yzi 

^liz) := eM-as,z-'l))x{z-'fiZ^T- 
Proposition 4.7. The following equations hold true: 

Li*i(z) = z^i{z) ^i{z)Ri = z^i{z) 
L'^^Uz) = z-'^*2iz) $2(^)^2 = z-^^*2iz). 



(23) 

(24) 

(25) 
(26) 



(27) 
(28) 



Proof We will just prove the last equation, all the other ones are proved in a similar way. 
Prom the very definition we have 



mz)Rj = z-'mUz) ^ [x{z-'r{Z,'VR^ = z-'[x{z-')nz^') 
[xiz-')fA-\Z^r ■ 



-^2 — ^2^ ^{Z2 
The proof of the following proposition is straightforward 



1\T 



R2 — Z2 AZ2 



□ 
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Proposition 4.8. The Lax equations h21\) and /i22\) are the compatibility conditions of the 
eigenvalue equations \21^ and Ii2^) with the following equations: 

(29) 
(30) 
(31) 
(32) 



Remark 4.9. Actually our Lax equations as well as equations (f27|l - (f32]l can be deduced from 
the equations of multicomponent 2D-Toda |10j . 



5 Recursion relations for matrix biorthogonal polynomials on 
the unit circle. 



In order to generalize the scalar theory we have to construct an analogue of the recursion 
relation given by Proposition 13.11 Recursion relations for matrix orthogonal polynomial on 
the unit circle are already known, see [TH] and [22j . Here we slightly generalize to the case of 
matrix biorthogonal polynomials on the unit circle. 
We define the following important n x n matrices: 



Definition 5.1. 



K 



N 



SC{Ti 



N+lJ 



SC(T^ 1 1) 



where SC denote the nxn Schur complement of a block matrix with respect to the upper 
left block (see for instance [21J). For example we have 



SCiTj 



7 



(0) 



7 



V^(-i)y 



(here and below Tj^'^ := (T^) ^ and similarly for /i^, /i^ and T^). 
Proposition 5.2. Monic biorthogonal polynomials such that 



^ p^r p(l)r . 



p(l)l p(2)/ , _ . , « 



are given by the following formulas: 

>(l)r 



P 



N 



SC 



7 



(N-l) 



\ I 



7 



i-N+l) 



Z 



,(-1) 
zH / 



/ 7 



(0) 



7 



(-A^+l) 



SC 



z'H J 



13 



p{i)iN ^ sc 



/ ^(0) 



( 7^°) 



V I z\ 



1 



(AT-l) J X 



7(0) z^-^I 

^(iV-1) ^(1) X 



^(0) ^(iV-1) 



Proof We will just prove the first formula, the second one is proved similarly. 
First of all we have VO<m<A^ — 1 



In the same way VO<m<-/V — 1 

f(^z-^i-h^^^ 7«)r- 



\z-^+HJ 



1 7(^)2; 



Finally 



< pS^^{z),pS^^{z) >r=< Z^1,P^^^^ >r 



h 



N- 



This completes the proof of the first formula, the second one is proved similarly. 



Remark 5.3. Note that, imposing 



we obtain biorthonormal polynomials. 
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Now we will write a long list of relations among this polynomials and reflection coefficients. 
Given any matrix polynomial Q{z) of degree n we define the associated reversed polynomial 
as 

Q(z) =z"Q*(z). 
The reflection coefficients are defined as follows: 

Proposition 5.4. The following formulas hold true: 

pWl _ ^p(l)^ _ p(2)r- /oo^ 

5{2)r 5(2)r _ r p{l)l (oa) 

pWr ^p(l)r _ p(2)l r /o.N 

p{2)l p(2)l _ ^p(l)r I .of5^ 

^N+l = ^-^At'' (I ~ yTV+l^A^+l) + (37) 
-fAT+l = 2^(1 — x5v+iyAr+l)-^Af + ^A^+l-^Af+l l*^^) 
-Pat+I = (I ~ 2/7V+l^Af+l)-^Ar'* +yAf+l-fAr+l (^9) 

-Pat+i = -^jv^ (I ~ ^iv+i^JV+i) + -fjv+iyjv+i (40) 

•^TV^A'" ~ ^TV-^TV (41) 

y]v/i5v = /^?vy^ (42) 
^N^^N+i = I ~ ykf+i^JV+i (43) 
^Wi^jv' =I-^Wiy7V+i- (44) 

Proof The first four formulas are proved observing, for instance for the first case, that 
V 1 < i < we have 

=< pl^f, - zpW', z^i >,=< ^ zn >i 

so that Pj^^^'-^ — zP^^"*' and P^'*'" must be proportional. Setting z = you also find the constant 
of proportionality. In particular, when proving (j34p and (j36p . we find a formula and then we 
have to take the reversed one. (f37]) is proved substituting (f36]) into (f35]) and similarly for 

(IMl),®,®. 

(|4ip and (j42p are proven respectively observing that we have 



and 



^ 6(1)' pClK . _ , p{l)l p(l)r . 
^ Af ' A'" — ^ A" ' Af 



p(2)r p(2)« _ p(2)r p(2)i 
A" ' A" ' — A" ' Af ' 



and then doing explicit computations. Finally (|43p is obtained rewriting (j37p as 

pWr p(l)r 

Ijv+i _ fA^.j ,/ t'' ^-urp^^^hv^ 
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multiplying from the left for -P^'^'t and then taking the residue. (j44p is proved similarly. □ 

Now we define two sets of block matrices {jC^}Ar>o and {C''j^}n>o- They will have, in the 
matrix case, the same role played by {Cn}n>o in the scalar case. 



Definition 5.5. 



N 



f -^WA . (46) 



Corollary 5.6. The following block matrices recursion relations are satisfied 

Proof These are nothing but ([33]), ([MI), ([35]) and ([36]) • □ 



6 Explicit expressions for Lax operators and related semidis- 
crete zero-curvature equations. 

Using our recursion relations we want to find explicit expressions for Lj and Ri in terms of our 
reflection coefficients x^j^,x^,y'j^,y^. First of all we underline a remarkable symmetry that will 
allow us to reduce the amount of our computations. Doing the following three substitutions 

z 1-^ 
t 1-^ —s 

s ^ —t 

we obtain immediately the following proposition. 

Proposition 6.1. Under the symmetry above the dressings, the orthogonal polynomials, the 
Lax operators and the reflection coefficients change as follows: 

rpr ^J,l rpl ^ j,r 

Si 1-^ S2 ^ Z^^ 

Li ^ R2 L2 ^ Ri 

^ (Q(2)r)* Q{2)1 ^ ^Q(l)ry 
p(l)« ^ (p(2)r)* p{2)l ^ (p(2)ry 

y^ y'' 
hi ^ hi 
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Hence we can write just the left theory and we will have the right one as well. Actually 
every computation made above for the right theory can be deduced from the left theory and 
this symmetry which will be called in the sequel t — s symmetry. In the theorem below the 
symbol rijf=Ar+2 means that the terms in the product must be written from the smallest index 
to the biggest, going from left to right. The symbol nj=1v+2 means that the product must 
be taken in the opposite direction. 

Theorem 6.2 (Lax operators for block Toeplitz lattice). Lax operators Li and Ri are ex- 
pressed in terms of reflection coefficients according to the following formulas: 
M N > M>-1 

M- 

{Li)nm+i = -xn+i{ n (i-yj4))yM+i (49) 

j=N+2 

M- 

{R2)n,M+1 = -yN+l{ n (50) 
j=N+2 
AI+ 

(L2)Af+l,7V = -^M+l4m( n i^-yj^'j))yN+lhN (51) 

j=N+2 

M+ 

(i?i)M+i,JV = -/iM+iyi/+i( n (l-^iyj))^iv+i^jv (52) 



j=N+2 



Moreover 



{Li)n,n+i = {R2)n,n+i = I (53) 
iL2)N+i,N = h^N+i^^ {Ri)n+i,n = h'^N+ih^' ■ (54) 

Proof (j53p and (j54p follow trivially from the expressions of dressings Si, Zi. Infact, because 
of our normalization, we have that 5i and Z2 are equal to the identity matrix plus a strictly 
lower triangular matrix which is what stated in (j53p . (j54p is obtained observing that the block 
on the main diagonal of 5*2 and Zi are given respectively by matrices /i^ and h^. as follows 
from Proposition 14.21 and Remark 15.31 Now let's begin with ()49p and (jSOp : the important 
point is that we have 

^1 =exp(^(t,zl))p(^)'. 
Hence as done in [T3] we can find that V > Af > — 1 we have 

(-^-i)Ar,M+i = -x^n+i^'n^JI+iVm+i- (55) 

Infact ViV > Af > -1 

^ p(l)' _ ^p(l)' p(2)« _ ^p(2)« ^ _ _ ^ ^p(l)' p(2)« . _ 
- < Pn+1 + • • • + {Li)n,M+iPm+1 + • • • ) Pm+1 >l= - {Li)n,R4 +1^^1+1 ■ 

On the other hand, using recursion relations, I also have VA^ > Af > — 1 

^ p(l)« ^pWl p(2)« ^p{2)l I p(2)r . I .Tp{l)r 

^^A^+l y J> Z K^N ) 1\^)^M ^~ I Vm+I 



J ^ pi'2)rN-Mpil)r I _ I ,r / 
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and comparing them we find (|55p . Now we use t — s symmetry to simplify this expression. 
We obtain 

iR2)N,M+l = -Un+I^N^M+I^M+I- 

Now we apply several times the recursion (ji^ to the piece ^jv^m+i S^t (f50]l . (H9]l is ob- 

tained using t— s symmetry. For ()5ip and (I52p we start defining i?i such that zP^^^"^ = P^^^'^Ri; 
then we will have Ri = h'^Rih~'^ and computations for Ri is carried on similarly as for Li. □ 



Remark 6.3. Our equations (P^ . (fSn]) . ([^T]) and ([5^ extend to the matrix biorthogonal setting 
the equations written in [22j for matrix orthogonal polynomials (see equation (4. 2), (4. 3)). In 
that article properties of M are applied to study some problems in computational mathema- 
tics (multivariate time series analysis and multichannel signal processing) and no relation is 
established with the Lax theory and integrable systems. 

The theorem above describes the block-analogue of Toeplitz lattice; we are now in the 
position to prove the analogue of Theorem 13.41 

Theorem 6.4. Block Toeplitz lattice flow can be written in the form 

for some block matrices M.^^ ^, M.^^^ ^,M.\. ^ depending on the matrices {x^, yj, x^, yj} 

and the spectral parameter z. 

Proof As we did for the scalar case we prove it just for t times. The relevant equations 
linking biorthogonal polynomials with wave vectors are 

^i(z) =exp(e(t,zI))P«'(z) 

$i(z) =exp(e(t,zl))p(i)'^/i-" 
= exp{-as, z-H)){h-'fdi[z-'])iP^^^Y. 
Trivial computations give the following time evolution: 

dt„P^'> = -d{[z]){R-,)^d{[z-'])p('> 

dt^p(^y = p^^>{h-''{R^).h' + h'^'idt^h"-) - z"i) 

the last two can be simplified giving 

5t„p(i)' = (L^)+P(i)' - z"ip(i)' (58) 
dt^P'^'^^' = -d{[z]){R^),di[z~'])P^^> (59) 

= p(i)'^(/i-^'(^5^)— /i'' - ^"i) (60) 

^ -pi2)i(^d{[z-'])h-'{u{)++hU{[z])y (61) 
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where (ii") means the lower triangular part including the main diagonal and (-L")++ means 

the strictly upper diagonal part. This simplification can be obtained evaluating the terms 
{dt^h"^) and h~\dt^h!')) through Sato's equations. Also they can be obtained observing 

fl)r (2)/ 

that ^'^d Pjy are monic so that the derivative of the leading term is equal to 0. Then 

the proof is obtained as we did in the scalar case using forward and backward recursion 
relations ([Ml), dSS]), ([391) and (00]). □ 

Corollary 6.5 (Non-abelian AL semidiscrete zero-curvature equations). The matrices 
and C}^ satisfy the following time evolution 

du/sA = K/s.,n+An - 4-^W.„n (62) 

= /:;a^[^/,,,„+i - (63) 

Proof These equations are nothing but compatibility conditions of recursion relations (j47p 
and (HE]) with time evolution ([56|) and ([57]). □ 



Remark 6.6. It should be noticed that, with respect to the equations originally written in [3], 
here we have two coupled non-abelian Ablowitz-Ladik equations. 

Example 6.7 (The first flows; non-abelian analogue of discrete nonlinear Schrodinger) . As we 
did for the scalar case we will compute the first matrices Ai^^^^^ ^ and use them to construct 

the non-abelian version of discrete nonlinear Schrodinger. We start with 7V4j ^; (I58p gives us 
immediately 

"ii^k — ^k+l ~ ^k+lVk^k ~ ^^k — 
^^k + ^k+l^k ~ ^k+lVk^k ~ ^^k — ~^k+iyk^k + ^k+l^k 

while we obtain immediately from (j59|) that 

0^*1 -^fc - zrLk'^k-i^k-i- 



Then we use recursion relation (I39p combined with 

/ifc/ifcli = (I - yUk) 

(this one comes from recursion relation ()4ip combined with t — s symmetry) to arrive to 



These computations give us 



I r I ■ 
'^k+lVk ^k+1 



A^Lfc =1 1 • (64) 

zyl -zl^ 



Exploiting t — s symmetry we can write immediately 



M's.^k =\ I • (65) 

-vl+i yl+iA , 
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Analogue computations for Ail^ gives us 



and 



a p(^y - pWr _ p{l}r I r _ pl{r) _ 
^h^k - Vk^k+l ^^k - 

^-f^fc "T ^k ~ Vk-^k+1 ~ ^-f^fc — ^k ~ ^k Vk-^k+1 



(here we started from ([HU]) and (fnT|) and we used recursion relations ([57j) . ([in]) . (jl^ . the last 
one combined with t — s symmetry). Then we arrive to 

(-yixl+i zyi' 
Ml,, = I (66) 

and using again t — s symmetry we also get 



^k+l 



^ -Vk+i 

M:,,, = 1 I . (67) 

-z -^kyk+i/ 

As we did for the scalar case we introduce times to and sq that give matrices 

,r/l _ fO 



Then we construct the matrices 



Z -^k+ltlk -^k+l ^ -^k 

zyl - yl+i yl+ixi + l-zlj 

'z-H - I - y^x^+i zy\ - y{^^ 

4+1 - ^"^4 4^1+1 -z\ + \, 
associate to the time r = ti + si — to — so- Semidiscrete zero-curvature equations 
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are equivalent to the systems 



( 1 1 0/1/ / r J 1 r I 



drxi 



Vk+i + 2yfc - vl-i + vl+i^ivl + yl^'kvl-i 
2x1 + xi 



X 



k+l 



^k-l 



ry* ^ ■? / ^ /y»^ ry*^ 11^ ^y* ^ 

•^k-lilk-^k -^kyk-^k+l 



(70) 



(71) 



Note that both of them are equivalent to the discrete matrix NLS as written, for instance, in 
[6] . Using (I70p and ()7ip together we perform the reduction to the hermitian case in a different 
way from [B]. First of all we rescale r i— > ir and then we impose 



r 

Vk 

yi 



±(4)* 
±(4)* 



Note that this reduction (with the sign plus) corresponds to studying the theory of matrix 
orthogonal polynomials on the unit circle as described in [19j and [22], hence it is very natural. 
This reduction gives us the two coupled equations 



-idrxi 



■^k+1 '^■^k ^ -^k-l ~r ^k+l\^k) ■^k + -^k\-^k) -^k-l 



(72) 



lUrXf^ — JLf^_^_^ ZX^ -t- X^_^ -|- X^_-^|^X^; Xf^ -\- Xf^^Xj^j X^_,_j^ 

already studied in ^ and generalized in 0. 

Remark 6.8. In [16] the authors studied finite gap solutions of the Ablowitz-Ladik hierarchy. 



It could be interesting to generalize their results to the non-abelian case. We will consider 
this problem in a subsequent publication. 
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